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Abstract 

We extend our previous work devoted to the computation of the next-to-next-to-leading order 
spin-orbit correction (corresponding to 3.5PN order) in the equations of motion of spinning compact 
binaries, by: (i) Deriving the corresponding spin-orbit terms in the evolution equations for the spins, 
the conserved integrals of the motion and the metric regularized at the location of the particles 
(obtaining also the metric all-over the near zone but with some lower precision) ; (ii) Performing the 
orbital reduction of the precession equations, near-zone metric and conserved integrals to the center- 
of-mass frame and then further assuming quasi-circular orbits (neglecting gravitational radiation 
reaction). The results are systematically expressed in terms of the spin variables with conserved 
Euclidean norm instead of the original antisymmetric spin tensors of the pole-dipole formalism. 
This work paves the way to the future computation of the next-to-next-to-leading order spin-orbit 
terms in the gravitational- wave phasing of spinning compact binaries. 
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I. INTRODUCTION 



Coalescing binaries of compact objects — neutron stars and/or black holes — are one of 
the most promising sources for the detection of gravitational waves (GW) with the advanced 
versions of the ground-based interferometers LIGO [l|, 0] and VIRGO js], 0], with GEO-HF 
[sl which will start taking data in the next years, and with the cryogenic detector KAGRA 
joining the network in a relatively near future |6|. Further ahead, LISA- type space-based 
interferometers 0, Isf, which will significantly increase the accessible region of parameter 
space, will be able to detect supermassive black-hole binaries with a very high signal-to- 
noise ratio. Successfully extracting the possibly very weak signal from the noise or estimating 
the parameters of the source with good accuracy and precision can be achieved by using 
matched filtering techniques, provided that the waveform is priorly correctly modeled (see 
e.g. Refs. [q, 10|). The post-Newtonian (PN) approximation scheme enables the computation 



of precise / accurate templates for the inspiralling phase of the coalescence of compact binaries 
ll|. 

For non-spinning objects, both the dynamics of the system and the waveform phase have 
been derived up to the 3.5PN order [l2|-[l3| Fl In the last years, an important effort (motivated 



by astrophysical observations |l5l-ll9l|) has been undertaken to extend these results to the 



spinning case, up to the same accuracy. In the present work, we investigate the higher-order 
binary dynamics focusing on the spin- orbit effects, which are linear in the spin parameter and 
numerically the most important to be taken into account at leading order. Here by spin, we 
mean the intrinsic (classical) angular momentum 5* of the individual compact body, rescaled 
in our convention by a factor c as 5* = c^true = G'm^^jyX, where x is the dimensionless spin 
parameter equal to 1 for maximally spinning objects, so that S appears to be formally of 
Newtonian order and all powers of 1/c are kept explicitly. Adopting this power counting, 
the leading order spin-orbit and spin-spin contributions to the dynamics appear at 1.5PN 



and 2PN respectively [2Ch25| . while the next-to-leading corrections are of 2.5PN [26H31| and 



3PN order [32N34j respectively. Note that the spin-spin couplings between different spins (1 
and 2) are actually known at 4PN order 35, 3y. 

In a recent work [s^] B we computed the next-to-next-to leading order spin-orbit correction 
at order 3.5PN in the equations of motion. We used a direct post-Newtonian iteration of the 
Einstein field equations in harmonic coordinates and proved the equivalence of our result with 



the one obtained in Ref. [38| using a Hamiltonian approach in ADM coordinates. We have 
further computed the associated conserved energy, verified the manifest Lorentz invariance 
of the equations of motion, and found agreement in the test mass limit with the motion of 
a spinless test particle around a Kerr black hole, as well as that of a spinning test particle 
around a Schwarzschild black hole. Our calculation was based on the description of the two 
compact objects as spinning point particles within the framework of the pole-dipole effective 
model. We refer to Section II of Paper I and references therein for a detailed presentation of 
the formalism on which the present work also relies. Our purpose here, building on Paper I, 
is threefold: 

1. Deriving the next-to- next-to-leading order spin-orbit corrections to the evolution equa- 
tions for the spins, the conserved integrals of motion and the metric regularized at the 
location of the particles. We also obtain an expression for the spin-orbit part of the 

^ As usual the nPN order includes corrections up to the relative order in a power expansion in the 

inverse of the speed of light. 
^ Hereafter this work will be referred to as Paper I. 
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metric valid all over the near-zone up to the highest available order. Knowing explic- 
itly the spin contributions in the near-zone metric is important for applications such 
as the numerical study of the accretion disk or jet dynamics around spinning back hole 



binaries [39|, and the on-going comparison between PN predictions and the numerical 



calculation of the self force acting on a particle orbiting a spinning black hole [40|, |41 



Expressing the results in terms of spin variables with conserved Euclidean norm instead 
of the original antisymmetric tensors of the pole-dipole formalism as used in Paper I. 
This choice simplifies the dynamics since the spin evolution equations take the form of 
ordinary precession equations, and is motivated by our final goal, namely to compute 
the next-to-next-to-leading spin-orbit effects in the GW orbital phase. The evolution 
of the GW phase will be indeed obtained via an energy balance argument requiring 
specifically the use of the conserved-norm spins. 



3. Performing the reduction of our results to the center-of-mass frame and to the case 
of quasi-circular orbits. This is of particular interest for data analysis purposes since, 
generally, orbits will have been circularized by emission of gravitational radiation by 
the time the binary enters the observational band. Notably, we obtain the 3.5PN 
expression including spin-orbit terms of the binary's energy as a function of the orbital 
frequency, which is one of the main ingredients entering the construction of the various 
approximants used for the evolution of the GW phase. 

The paper is organized as follows. In Section |lTl we present the construction of the 
conserved-norm spin variable that we use throughout the paper as well as the associated 
precession equation including the next-to-next-to leading spin-orbit effects, written in a 
general frame. The other general-frame expressions being too long to be shown here, we 
directly move on to the presentation of our results in the center-of-mass frame in Section 
mil The reduction to quasi-circular orbits is performed in Section IIV[ and the Section |V] 
contains our conclusions. In Appendix [X] we compute the correspondence between our spin 
variables and those of Ref. jlsf, whereas Appendix [B] gives the explicit relation between our 
spin variable and the antisymmetric spin tensor of Paper I. Appendix [C] gives the general 
frame expression of the center-of-mass position at next-to-next-to leading spin-orbit order. 
We relagate in Appendix [D] the lengthy expression of the total angular momentum. In 
Appendix [E] we show that spin-orbit terms at 3PN order are pure gauge. 

Most of our computations were achieved by means of the package xAct, which handles 
symbolic tensor calculus within the scientific software Mathematica® (42| . 



II. PRECESSION EQUATIONS OF COMPACT BINARIES 
A. Definition of the constant magnitude spin 

While in Paper I we found most convenient to work with the space components of the 
spin tensor, in the present article we shall introduce new spin variables, denoted Sa (with 
a = 1, 2, 3), which are designed to have a conserved Euclidean norm. Using conserved-norm 
spin vector variables is indeed the most natural choice when considering the dynamics re- 
duced to the center-of-mass frame or to circular orbits. Their evolution equations reduce, 
by construction, to ordinary precession equations, and these variables are important when 
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studying the gravitational waves emitted by a quasi-circular binary because they are secu- 



larly constant [43|, |44 



For convenience, we change our notations for the spin variables with respect to Paper I. 
For the spin tensor variable of Paper I we now use the notation: 

^ — '-'mbfb • \^-^) 

Let us recall from Section II in Paper I that this variable satisfies the covariant spin supple- 
mentary condition (SSC), namely S^^p^ = where denotes the linear four-momentum 
of the particle. The covariant spin vector 5^ associated with the spin tensor is therefore 
defined bjH 

S""" = ^ e^^P^^S^ , (2.2) 

y/—g m 

where e'^"'"^ denotes the Levi-Civita symbol; we know that both p^p'^ = — and S^S^ = 
are conserved along the trajectory: m = const and s = const. Working at linear order in the 
spins, the linear momentum agrees with the normalized four velocity, p^ = mu^ + 0{S'^), 
and we simply get that: 

^ = 0{S') , (2.3) 
dr 

where D/dr = u'^'S/^. Thus the spin covector is parallel transported at linear order in spin. 
We can also impose that the spin should be purely spatial for the comoving observer: 

SX = OiS^) . (2.4) 

From now on, we shall omit writing the (9(5'^) remainders. 

A standard, general procedure to define a constant-norm spin vector consists in projecting 
onto some orthonormal tetrad e^^, i.e. a tetrad that satisfies g^iv^a^^p = Va/s, which 

leads to the four scalar components 

Sa = e/S^ . (2.5) 

If we choose for the time-like tetrad vector the four-velocity itself, Cq^ = m^, the time 
component tetrad projection 5*0 vanishes because of the previous orthogonality condition 
S^u^ = 0, and the spatial components Sa (with a = 1,2,3) define a constant-norm spin 
vector. Indeed we have seen that S^^S^ = = const is conserved along the trajectory. 
Because S^u^ = we can rewrite this as '^^'^S^S,, = s^, in which we have introduced the 
projector 7^^^ = g'^'^ + u'^u'^ onto the spatial hypersurface orthogonal to u^^. Now, from the 
orthonormality of the tetrad and from our choice Cg^ = u^, we have 

= ^"'e/e,^ . (2.6) 

Therefore the conservation law ^^'^S^Sy = becomes 

d'^'SaS, = , (2.7) 

which is the relation defining a Euclidean constant-magnitude spin variable Sa- 



For convenience in this Subsection we pose c = 1. 
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However, the choice of the spin variable Sa is still somewhat arbitrary, since a rotation 
of the tetrad vectors can always be performed. Here, in order to fix it, we shall achieve 
a construction equivalent to that of Ref. [29[, which presents the advantage of providing a 
unique determination of these redefined variables in a given gauge. This results in a definition 



which differs from the one adopted at the previous order in Refs. [28|, |4J]. In Appendix R] 
we shall study the correspondence between our constant magnitude spin variable and the 
ADM one, used in Refs. Q Hi. 



To uniquely fix the tetrad vectors we proceed in the following way. Consider the spatial 
covariant components of the projector tensor 7^1,, namely '-fij = Qij + UiUj, which as we see 
from Eq. (12. 6p can also be written in terms of space tetrad vectors as = 5°''^ eaiCbj ■, where 
the "mixed-components" tetrad vectors read Cai = e^^'-^i^. The problem is thus to find a 
prescription for defining the 3x3 matrix Cai starting from the 3x3 matrix 7jj. Clearly, if 
Cat is symmetric in its two indices, namely Cai = (in which we are actually exchanging a 
covariant spatial index i with the Lorentz spatial index a), then the matrix Cai can be viewed 



as the square root of the matrix 7jj. Now, as pointed out in Ref. 29|, a positive-definite 
symmetric matrix such as 'jij has a unique symmetric positive-definite square root Cai, which 
therefore satisfies 7jj = S'^^Ciaebj- Thanks to this lemm£0 we have a mean to define the tetrad 
in a "canonical" way, by adopting for Cai the unique symmetric positive-definite square root 
of 'jij. We find 

e/ = [Y' - Y'v^) ea^ , (2.8) 

where we have used Cao = —v^Cai with = /vP denoting the ordinary coordinate velocity. 
This, together with = u^-, defines completely and uniquely the tetrad and therefore the 
constant-magnitude spin variable 

Sa = e/S, = {e: - v^e^') . (2.9) 

The evolution equation (12. 3p for the original spin variable S*^ now translates into an 
ordinary precession equation for the tetrad components (12.91) . namely 

^ = ^aSb, (2.10) 

where the precession tensor is given in terms of the (covariant) derivatives of the tetrad by 

Be'' 

n: = -e,^^. (2.11) 

The antisymmetric character of the matrix Qj', which is just made out of the usual Ricci 
rotation coefficients u°'J' = —e°'^'WuG^^ = — w''^", namely = v'^u:°'J' where v'^ = /vP, 
guaranties that Sa satisfies an ordinary precession-type equation, i.e. 

A G 

— = X S , (2.12) 

dt ^ ^ 

where we denote S = (Sq), Q, = (Qa) and pose Qa = —^^abc^^'^- As a consequence of (I2.12p 
the spin has a conserved Euclidean norm, = [cf Eq. (12. 7p . 



^ The existence of a symmetric positive-definite square root of a symmetric positive-definite matrix comes 
from the fact that the latter can be diagonahzed and has only positive eigenvalues. The unicity is 
easily proven after remarking that the original matrix and its square root commute and therefore can be 
simultaneously diagonahzed. 
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The latter construction is completely equivalent to the one of Ref. H. Indeed the 



orthonormality condition g^^^S^Sy = is rewritten in |29| as G^^SiSj = s^, with an effective 
metric G*-' = — 2g^^'^v^^ + g'^^v^v^ , and the conserved norm spin variable is defined there 
by S** = H^^Sj, where H^^ is the unique symmetric and positive definite square root of 
the matrix G^^ . Replacing the metric g'^'^ by the projector 7^^*^ = g^^'^ + u'^u'^, and using 



7' 



fJ,U 



6 eJ^Cfj", we obtain the correspondence with our previous tetrad formalism: 



(2.13) 



Notice that the effective metric G*-' is actually the inverse of the spatial projector = 
dij+UiUj, i.e. we have G^^jjk = SI. From fl2.13p we get immediately the correct factorization 



Indeed, the matrix such that G^^ 



^kijj ijj 3 ^Yie notation of 29 , can be written as 



V e„ 



(2.14) 



If we replace the tetrad components by (12. 8p in this expression, we find that 
Since the symmetric matrices Cai and G^^ are, respectively, the square root and the inverse 
of the matrix 7jj, we obtain that the matrix H^^ is actually the inverse of the matrix Cai, 
namely Hjchj = Sab, which shows that it is symmetric and positive definite, as required. As 
a result, we find that the definition of Ref. 29 



(2.15) 



is in fact the same as our definition (12. 9p . Both procedures explained above are therefore 
equivalent, but the present approach based on a tetrad is more insightful. 



B. Evolution equation for the conserved spin vector 

From now on we shall stick to our definition (12. 9 p of the constant magnitude spin, and 
shall no longer make any distinction between the spatial indices a,b, - ■ ■ , which were intro- 
duced above as tetrad indices, and the ordinary spatial indices i,j,k,--- = 1,2,3, which 
will be raised and lowered by the Euclidean metric. The 3.5PN gravitational field will be 
decomposed for convenience in terms of eight metric potentials denoted by V, Vi, Wij, Ri, 
X, Zij, Yi and T, as introduced in Eqs. (3.1) of Paper I: 



900 



X + ViV, + 



32 



9oi 
9ij 



16 



48 
1 



Sij 



1 1 



16 



6 



+ ^[T- ^-VX + R,Vi - ^-VV.Vi - _ ) + O ( -L 



04 



(2.16a) 
(2.16b) 



(2.16c) 
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Each of these potentials is a retarded solution of a fiat-space wave equation sourced by 
components of the stress-energy tensor and appropriate lower order potentials. We refer to 
Eqs. (3.4) of Paper I for the explicit expression of these source terms. 

Our first result is the explicit expression of the precession tensor in terms of the 
metric at the 3PN relative order. For a single spinning particle described in the pole-dipole 
approximation and moving in a fixed background, the evolution equation for the spin is 
given by Eq. (12. Sp . The metric and associated Christoffel symbols therein are replaced up 
to 3PN order by the potentials V , ■ ■ ■ , Yi, and the spin variable is changed to conform with 
the definition (12. 9p so that we can read off the precession tensor from Eq. (I2.10p . In practice 
we use the following explicit expression in terms of the matrix defined by Eq. (I2.14p : 



dt 

where the matrix K^^ is a combination of velocities and Christoffel symbols: 

,k 



(2.17) 



Oi^ ^ '- ik 



(2.18) 



We then obtain the anti-symmetric precession tensor up to 3PN order as 
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-Wi^v^'d^y^^ + m'^'kd^y^^ - v^v^v^'diW^u + v^v^'dtW^^^ ■ (2.19) 

This quantity is to be evaluated at the location of the considered particle, at which point 
all the potentials V , Vi, ■ ■ ■ , that parametrize the metric are diverging. For a system of 
particles, the metric is the one generated by the stress-energy tensor of the system of the N 
particles itself. We consider the case N = 2 and compute the precession tensor (I2.19P at the 
location of, say, particle 1, the velocity appearing in (I2.19P being thus the ordinary velocity 
of particle 1. 

The evaluation of (I2.19P at 1, i.e. at point yi, is made meaningful through Hadamard's 
regularization 45|, consistently with our computation of the 3.5PN spin-orbit acceleration 
in Paper I. It must be understood in Eq. (12.171) that each term is regularized before taking 
the products and the time derivative. Nonetheless, although Hadamard's regularization is 
"non-distributive" in the sense that {FG)i ^ {F)i{G)i in general, we checked that, at the 
order considered here, treating it as "distributive" makes no difference. Note that the time 
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derivative operation does not commute with the regularization operation at 1, and we have 
generically for singular functions F in the class considered in Ref. Bill 



^(F)i = (9,F)i + K9,F)i, (2.20) 

where (G)i represents the value of G at particle 1 position in the sense of the Hadamard partie 
finie. In order to present a closed-form expression for in terms of the metric potentials, 
we first applied the total time derivative there according to the Leibniz rule on individual 
monomials composing , applying the distributivity ansatz [i.e. {FG)i = (F)i(G)i] for the 
products. We next replaced the accelerations by their expressions in terms of the potentials. 
For the time derivatives of quantities regularized at 1, we resorted to Eq. (12.20 p . Finally, the 
partial time derivatives of the potentials were eliminated in turn by means of the identities 



(3.28) of Ref. j46|, which are equivalent to the harmonic gauge condition. 

Since we are working at linear order in the spins, only the non-spin parts of the metric 
potentials enter the computation of the matrix . Most of those contributions are the 
same as those required for the 2PN equations of motion without spinjl There are only two 
genuine 3PN potentials: One of them, Zij at Newtonian order, has the same structure as 
Wij] The other one, Yi, which enters the term —IGd^'^Y^^ in Eq. fl2.19p . shows a higher order 
of non-linearity (in powers of G). Only its regularized value can be computed, using dimen- 
sional regularization in principle, as was done for the 3PN equations of motion without spin 
obtained in 0. Like for th; term S^^id.^Y,) 1 appearing in the equations of motion (see 
Section V of Paper I), we find that the corrections coming from the dimensional regulariza- 
tion exactly cancel out because of the antisymmetrization due to the contraction with the 
spin tensor. Thus, like in Paper I, Hadamard's regularization is sufficient for our purpose 
here. The remaining 3PN metric potential, T, does not contribute. 

Due to the length of the expression, we relegate to Appendix |B] the relation between 
the conserved spin vector and the spin tensor in terms of the orbital variables derived from 
Eqs. (12. 2p and (12. 9p . We conclude this Section by giving the explicit expression for the 
precession equation of the conserved spin 1: 

^ = f^ixSi. (2.21) 
dt 

The vector Qi may be expanded at 3PN order in the form: 



= lo}™ + Afif ^ + + o(\] . (2.22) 



Except for the spin tensor, we use the same notations for the orbital variables as in Paper I: 
[uv) denotes the scalar product u ■ v = u^v^ and w = u x v the cross product between u 
and V, whose components are given by = e^^^u^v^. At leading order, we have 



sir = 

' 12 



G 



3 

- ni2 X vi - 2 ni2 x V2 



(2.23) 



^ This equation states that, formaUy, the Hadamard regularization commutes with the operator v^dp,. 
^ The non-spin part of the acceleration has the form a' — F'^ — dQ^/dt with — P' — w*; sec Eqs. (3.5) 
and Eqs. (3.7) in Paper I. 
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2PN _ G (2)^0,1^^ ^ p)n',''ml + ^'^nl''mm2] 



12 



r 



12 



while the next-to-leading order correction is given by 
where 

^ 'ill = ^12 X Vi 

L ^ 

+ ni2 X V2 3(r;,i2f2)^ + '2{viV2) - 1v\ 



(2.24a) 



^-^[nx2V2f + ^t'? - (fi^^2) + vi 



Vi X V2 



\nx2V\) - -[ni2V2] 



(2)fil'l = _ X Vi ni2 X V2 , 

(')\ 2 1 5 

^^-Ti^' = -- ni2 X vi - ni2 X V2 . 



(2.24b) 



These results can either be directly computed from Eq. (12.19P truncated to the appropriate 
order or obtained by a mere translation of the spin evolution equation given in Section VI 
of Paper I using the relation (IBip of Appendix [Bl Using Eq. fl2.19p . we get the next-to- next- 
to-leading order correction, new with this paper, to the precession vector: 



^3PN ^^(3)^0,1^^ ^ [^^)nfml + (3)f2i'Vim2] 



12 



12 



^12 



(2.25a) 



where 



= ni2 X vi 



^(^12^2)"^ - ^(?^12^^2)^^^l + + ^(?^12t'2)^(^'l^'2) - 7^^l(^^1^^2) + V2 

lb lo lo z 4 



-^[ni2V2) V2 + -^1^2 - {VlV2)V2 



-3{ni2V2f{viV2) + 2{viV2)vl 



+ Vi X V2 
1 



+ ni2 X V2 

1 



15 



("-12^2) + Q{ni2V2) V2 - 



9 



?^12^^l)^^l + -(?^12^^2) + (ni2^^l)t^2 



3 11 3 5 

■-{ni2Vl){ni2V2y - -(?^12^2)^l - -(^12^l)(flf2) + - (^^12^ 2) (f 1^'2) - -(^^12^2)^2 

Z o Z Z Z 



(^)rii'^ = ni2 X vi 
71 



27 29 67 

-Y(?^i2^^i)^ + 27(ni2t;i)(ni2t;2) - 14(ni2t;2)^ + -^v^ - —{V1V2) 

109, ,2 101, , , ^ 73, ,2 37 2 
ni2 X V2 — (ni2t;i) —{ni2Vi){ni2V2) + —{ni2V2) - y^i 



35 
"2 



4 



Vi X V2 



71, X 55, ; 

-y (ni2^^i) + y (^12^-2) 



(3)0°'^ = ni2 X vi 

+ ni2 X V2 



1 3 3 3 1 

-(?^12^^l)^ - -(?^12^^l)(?^12^^2) + 2(7212^^2)^ + -vf - -iViV2) - -vj 

2 2 o o 2 



1 5 
-5(^12^2)^ - 2^1 ~ 2(^1^2) + 2^'2 



Vi X V2 



5 1 



9 



(3)^0,3 



69 
31 



-— ni2 X vi + -ni2 X V2 , 
o Z 



4 



ni2 X vi 



17 



ni2 X V2 



11 



ni2 X vi - Y ni2 X V2 . 



(2.25b) 



III. COMPACT BINARIES IN THE CENTER-OF-MASS FRAME 
A. Reduction to the center-of-mass frame 

In this Section, we present our results in the center-of-mass (CM) frame, defined by 
the nuUity of the center-of-mass position G. We remind that the latter vector is directly 
linked to the conserved integral of the motion K associated with the boost invariance of the 
conservative part of the dynamics, through the relation K = G — Pt, with P being the total 
linear momentum of the binary. We shall display G in a general frame at the 3.5PN order 
in Appendix O However, it turns out that only the 2.5PN order is needed to reduce the 
results to the CM frame, since 3.5PN corrections that would matter for the reduction of the 
Newtonian part can be seen to cancel out based on symmetry arguments, as detailed at the 
end of Appendix O When working in the CM frame, it is convenient to use the standard 
mass parameters m = rrii + m2, Sm = nii — m2 and u = mi m2/m? (such that < z/ < 1/4), 



and to introduce the same combinations of spin variables as in Refs. [28|, |4J|, namely 



S = Si + S2, (3.1a) 

S^mf^--^V (3.1b) 

\m2 nil/ 

Notice that since we are working with different definitions of Si and S2, our variables S and 
S also differ from those used in Refs. 28|, |4^ ; however, 4^ introduced in their Section VII 



conserved- norm spin variables, which turn out to be identical to ours at IPN order in the 
center-of-mass frame. Neither S nor S are of conserved norm and, therefore, their evolution 
equations cannot reduce to precession equations; they can straightforwardly be written in 
terms of the precession vectors ili and 0^2 as 

^= f!^!irii + — ^2') xs + z/(ri2-Oi)xE, (3.2a) 

dt V m m J 

dS f^\^_^ ^ ^^^2) X S + (ris - rii) X S . (3.2b) 

1- 



dt V m m 

Here, we shall simply give the CM-frame expression for f2 

The vector G including spin-orbit effects at 2.5PN order was computed in Ref. 28 
Translating it to our spin variables and imposing G = 0, we find the following relations 
between the positions yi, y2 in the CM frame and the relative position x = yi — y2 and 
velocity v = dx/dt = vi — V2 = Vi2 (as well as spin variables). Posing 

m2 

yi = — x-hz, (3.3a) 
m 
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y2 



mi 



-X + z 



m 



(3.3b) 



we obtain 
2c^ m 




Gm 



V ^ rv 6m / 3 3,4 
— SXV + — < iy]v^n 



m V8 2 



1 3 \ , ,2 /19 3 \ 2\ 7, ,■ 
'^^l^]^™> +V"8'+2^j^ ln--(nt;)v 



+ 



4-2'" 



1 \ 9 „ Gm ( 5m 

2 / r V m L 



3 1 
(n, S*, f ) n — - (m;) nxS — - Sxv 



+ (-1 + 4i/) (n?;) nxS-(2 + z/)Sxv 




(3.3c) 



Note that the exphcit relation between the center of mass velocities vi, V2 and the relative 
variables v, x, as well as the two spins, can be obtained either by taking the time derivatives 
of Eqs. (13.31) and using the evolution equations, or equivalently by imposing that the explicit 
expression of the total linear momentum P up to 2.5PN order satisfies P = 0, which must 
hold in the center of mass frame, at that approximation level. 

For each quantity of interest, we only present the (linear in) spin part, which already 
leads to somewhat lengthy expressions, and refer to the literature for the non-spin terms. 
However, one should keep in mind that the spin part of a quantity in the CM frame is 
not the reduction to the CM frame of the spin part of the same quantity expressed in a 
general frame, because of the crucial additional spin contributions due to the replacement 
of Eqs. (13.31) into the non-spin part. 

B. Precession equation 

In the center of mass frame, the expression for the 3PN precession vector Q,i defined by 
Eq. ^TM) reduces to 
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We obtain O2 from Oi simply by performing the substitution 5m — )■ —5m. 
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C. Relative acceleration 

The acceleration including spin-orbit effects up to 3.5PN order was obtained in Paper I. 
The reduction to the CM frame yields 
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where we have indicated the contributions of spin-spin terms (not considered in the present 
work); the non-spin terms can be found in Ref. ll|. In the rest of this paper, we use the 
notation [u,v,w) for the mixed product e^^^u^v^w^ . For the spin-orbit terms, we have at 
leading order the standard result j22|. 
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At next-to-leading order, the results of |28| may be rewritten, with our choice of spin vari- 
ables, as 



"^B2.5PN 

s 



Gm 



Gm 



(3.7a) 



where 

(2-5)bi = n X S 
+ n X S 
+ S X V 



45 , ^3 





45 \ 




(-1- 


1 ly 






2 ; 





^ 5m / 9 
+ (n, L, V) V — 7:: - 6z/ I (nv) 
m \ 2 



5m , ... 5m f 3 , „ \ / . o 

-Ihu {nvY H h 12z/ I {nv)v^ 

m m \ 2 



3 

'2 



-- - 15u 1 (nvf + Uuv"^ 



S X V 



5m 
m 



3 
2 



— I -- - 9i^ 1 (nvf + 7u—v^ 



5m 
m 



12 



+ (n, T,,v) n 



m m 



+ (n, ^, v) n [30z/(nt;)2 - 24i/w2] , (3.7b) 



(2-5)b2 = (n, S, t;) n— ( 24 + — z/ ) + (n, 5, n (44 + 33z/) + n x S (-28 - 29u) (nv) 
m \ 2 / 

/ 31 \ , , „ , N „ Sm ( 19 , 
+ n X S — -12 z/ {nv) + S x v (-24 - 19z/) + S x v — -12 v \ . 



(3.7c) 



Finally, at next-to-next-to-leading order, we get 
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D. Conserved integrals of the motion 

As already verified in Paper I, the dynamics associated with the next-to-next-to-leading 
spin orbit terms at 3.5PN order is purely conservative [i.e. there exists a conserved energy 
at that order) and, in the harmonic coordinate system we are using, is manifestly invariant 
under Poincare transformations. Therefore, it admits the usual set of Noetherian integrals of 
motion associated with the Poincare group (energy, linear momentum, angular momentum, 
center-of-mass integral) . 

These integrals can be computed in a general frame from the acceleration given in 
Eqs. (6.7)-(6.12) of Paper I (translated to our conserved-norm spin variable) on the one 
hand, and the 3PN precession equation given previously in Section [11 Bl on the other hand, 
by using the method of undetermined coefficients as was done for the energy in Paper I. Note 
that the precession equation is needed at 3PN order only for the total angular momentum 
while the 2PN order is sufficient for the other quantities. 

However, a technical difficulty of the method of undetermined coefficients is that it must 
deal with the fact that a given expression does not admit a unique writing due to the 
existence of dimensional identities relating the different vectors of the problem. For this 
reason, we have chosen here a more direct approach consisting in translating the expressions 
obtained in Ref. j38| from a reduced Hamiltonian formalism in ADM-type coordinates and 
then verifying that they are indeed conserved as a consequence of our equations of motion. 
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Note that in the Hamiltonian formahsm, the total hnear momentum and the total angular 
momentum are trivially expressed in terms of the canonical variables. The translation to our 
variables is performed by applying to the ADM results the contact transformation derived 
at the required order in Section VII D of Paper I. For the angular momentum, we also need 
the relation between our spin variables and the canonical spin, IPN order beyond that given 
in Paper I. The derivation of this 3PN correspondence between both sets of spin variables 
is presented in Appendix |X1 

The resulting expressions in a general frame happen to be very long, hence we directly 
display them reduced to the CM frame, in which, by definition, only the energy and the 
angular momentum do not vanish. The latter quantity being still relatively lengthy, it is 
relegated to Appendix [Dl The energy in the CM frame takes the form 
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where the non-spin terms can be found for instance in Ref. For the spin-orbit contri- 

butions, we get at leading order 
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which differs from [23|, |28| due to the deviation in the spin variables beyond leading order. 
For the next-to-next-to leading order, we arrive at 
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E. Near-zone metric at the location of a particle 



The post-Newtonian near-zone metric up to order 3.5PN is parametrized in terms of the 
elementary potentials according to Eqs. fl2.16p . The next-to-next-to-leading spin-orbit effects 
correspond to terms of order 0{l/c^) in Qqq, 0{1/c^) in g^i and (9(l/c^) in gij. However, 
because of the complicated non-linear structure of the source terms, analytical solutions to 
the appropriate order can only be obtained for some of these potentials using the technique 
described in details in Sections IIIC and D of Paper I. As a result, the spin-orbit contributions 
to the near- zone metric evaluated at an arbitrary field point in the near zone is computable 
in analytic closed form only up to 0{l/c^) for g^o, 0{l/c^) for g^i and 0{l/c^) for g^j. The 
corresponding expressions will be given in Section [III F[ By contrast, the value of all these 
potentials at the location of each of the two point particles, defined by means of the pure 
Hadamard-Schwartz regularization procedure, has been calculated explicitly by using the 
tools presented in Section HID of Paper I. This has enabled us to obtain an expression for 
the regularized metric at the particle positions to the desired orders. In addition to the 
regularized potentials already obtained in Paper I, we shall now need the spin parts of (^)i 
to next-to- leading order, as well as (T)i to leading order, which we evaluate in a similar 
fashion. We define 
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The non-spin parts can be found in Ref. [48|. The spin-orbit contributions at 3PN order, 
namely the terms ~ 1/c^ in s9oo and ~ 1/c'^ in s9oi^ are actually pure coordinate effects. 
These terms can be eliminated by means of a gauge transformation as verified in Appendix 
[El (extending Ref. [i^). After reduction to the center-of-mass frame, we are led to the 
following spin parts. 
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5. Metric components qqi 
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3. Metric components gij 
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F. Near-zone metric in the bulk 



Let us now turn to the spin part of the near-zone metric up to orders 0{l/c^) for qqq, 
0{l/c^) for and 0{l/c^) for gij. As mentioned above, the wave equations defining the 
potentials at play up to these orders can be analytically integrated. We refer to Section V 
of jisj for a thorough presentation of the relevant method. The non-spin parts of the metric 
up to the same orders can be found in Ref. j48| . 

The metric at a given arbitrary field point x is most conveniently expressed using the 
distances to the two particles ri = |x — yi| and r2 = |x — y2|, as well as the unit vectors 
111 = (x — yi)/^i and n2 = (x — y2)/i^2 pointing from the particles to the field point. In 
order to avoid any confusion, let us recall that, as in the rest of the present Section IIIIl 
r = ri2 and n = ni2 respectively denote the radial separation between the two bodies and 
the unit vector pointing from body 2 to body 1. 

One key ingredient in the computation of the relevant potentials is the function g = 
ln(ri + r2 + r), which satisfies Ag = l/(rir2), where it is understood that the Laplacian 
operator acts on functions of x. In the spin part of the metric presented below, g appears 
within terms of the form {diijd2kg), where dn and d2i denote the partial derivatives with 
respect to y\ and 2/2- explicit form for such terms is straightforwardly obtained by 
expanding these derivatives but we shall keep them factorized in an effort to reduce the 
length of our expressions. Finally, we find that using the spin variables Si and S2 instead of 
their combinations S and S leads to much more compact expressions and, since the metric 
has to be symmetric under the exchange 1 O 2, we can then further shorten our formulae by 
only explicitly giving the contribution containing Si. Note that in practice, the operation 
1^2 implies the replacements 6m — )■ —6m, v — j- — v and n — )■ — n in addition to those 
where the indices 1 and 2 explicitly appear. 
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2. Metric components Oi 
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3. Metric components ij 
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IV. COMPACT BINARIES IN QUASI-CIRCULAR ORBITS 



A. Reduction to circular orbits 

We now present our results in the case where the orbit is nearly circular (if one neglects 
the radiation- reaction effects), i.e. has a constant radius and is planar apart from small 
perturbations induced by the spins that cause the orbital plane to slowly precess. Most 
compact binaries will be in quasi-circular orbits when they are observed by the gravitational 
wave detectors. Assuming this sort of motion drastically simplifies the long expressions of 
the previous Section, which become immediately much shorter. We may thus include also, 
for the sake of completeness and convenience, all non-spinning terms up to 3PN order. The 
results that we display are then complete up to 3.5PN order and linear order in the spins 
(remember that we neglect spin-spin contributions) Q 

When comparing with the various expressions in the literature for the lower order spin- 
orbit corrections, we have to keep in mind that different definitions of the spin variables 



yield different numerical coefficients. In particular, as Ref. [28| and Section VI of Ref. [44 
did not make use of conserved-norm spins, the results there differ from the ones that we find 
here at next-to- leading order. By contrast, our energy, angular momentum and precession 



vector agree with those given in Section VII of [4^, where conserved- norm spin variables 
are introduced, because these spin variables, although constructed differently, turn out to 
agree with ours when evaluated in the CM frame. 



For the kinematics of quasi-circular orbits, we use the same conventions as Ref. [28 
which can be summarized as follows. We introduce an orthonormal triad {n. A, €}, such 
that n = x/r, £ = Ln/|Ln| where Ln = /ix x v denotes the Newtonian angular momentum, 
and A = € X n. All the relevant CM quantities have been defined in Section [IIIAI We project 
out the spins on the latter orthonormal basis, as S = ^nii -|- S\X + Se£ and similarly for 
S. By definition of quasi- circular orbits, the orbital separation decreases as r = 0(l/c^), 
entirely due to the gravitational radiation reaction. The spins contribute to this radiation 
reaction only at 1.5PN order beyond the dominant term 0{l/c^), i.e. at 4PN order j43| . 
which is neglected here. We define the orbital frequency u and the precession frequency w 
by the relations u = X.h and w = £.X, so that the triad evolution equations read n = a; A, 
A = —un + uj£ and £ = —zxjX. Then, posing = zun, we have £ = zu x £ which is the 
natural definition for the precession vector -cu of the orbital planejl The acceleration may 
be then written in our moving basis as 

'^ = -uj\n+(ruj + 2ruj)X + ai£ + o(\] , (4.1) 
dt V c / 

in which we have set = ruow and neglected a term f = 0{l/c^). The spin contributions 
in the term proportional to A are, by the same argument as before, of negligible 4PN order 
since the spin terms in uj come with a factor l/c" and can thus be regarded as constant 
at the considered order, see Eq. (14. 2 p below and the discussion thereafter. As a result, the 
spins can contribute to the acceleration only through u and the quantity at 3.5PN order. 



^ Note that in order to reduce the results of the previous Section to circular orbits, the contributions from 
the non-spinning terms up to 2PN order, wh ich are not presented there, must be taken into account. 



® Wc have vo — — Wproc in the notation of Ref. 
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B. Equations of motion and precession vector 



In the case of circular orbits, the relativistic generahzation of Kepler's law relating the 
orbital frequency uj to the body separation r, is most conveniently expressed in terms of 
the post-Newtonian parameter 7 = Gm/{r(?) = 0{l/c^). It is derived by reducing the CM 
relative acceleration given in Section IIII CI We find 
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(4.2) 



This expression is clearly gauge dependent since it involves the separation r defined in 
harmonic coordinates. Notably, the presence of a gauge constant Tq in the non-spin part 
at the 3PN order attests the occurrence of a pure gauge effect which disappears from the 
observables of the problem (see e.g. [lH). As we stated before, the spin part of u is constant 
at linear order in spin and up to 3.5PN order, for r is constant at the 2PN approximation 
and so are the spin components 5"/ (with A = 1,2): indeed, the precession vectors ^Ia are 
proportional to £, see Eq. (14. 5p . hence dS'//dt = —wS^, the right-hand side of the latter 
relation being quadratic in spin since w is proportional to the spin, see Eq. (14. 4p below. 

When computing physical quantities, it is useful to eliminate the coordinate distance r (or 
equivalently 7) by reexpressing it in terms of u or, more conveniently, of the post-Newtonian 
frequency-related parameter x = (Gmco/c^y^^. Inverting Eq. (14. 2p . we obtain 
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For the component of the acceleration along £, which vanishes in the absence of spins, after 
inserting the explicit value of v = rn + rcuA and eliminating r in terms of x by means of the 



24 



above formula, we find (recall that ai = ruw) 
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The last quantity that we need in order to complete the evolution equations for 
orbits is the precession vector fli. From Section [IIIBI we readily get 
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C. Conserved integrals of the motion 



Let us now provide the expressions of the energy E and of the orbital angular momentum 
L reduced to the case of circular motion. Applying to the CM expression (13. 9p the same 
reduction procedure as for the computation of in the previous Subsection yields 
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where the non-spin terms are taken from Ref. [llj. Notice that since is a physical observ- 
able, the gauge constant Tq has canceled out as expected. For the total angular momentum, 
we obtain from the CM expression provided in Appendix [Dl 
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The non-spin parts can be found for instance in Ref. 50|. At the next-to-leading order in 
spin-orbit effects, for both E and L, we recover the results of Section VII of Ref. 4J], which 
are given in terms of the conserved-norm spin variables defined there. 



D. Metric regularized at a particle location 

Finally, we give the expression of the spin part of the near-zone metric regularized at the 
location of body 1 for circular orbits: 
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V. CONCLUSIONS 

Building on previous work (Ref. jl^] , referred to as Paper I) we have presented complete 
results valid at the next-to-next-to-leading order, corresponding to 3.5PN order, for the spin- 
orbit contributions to the dynamics of spinning compact binaries. We notably obtained the 
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spin-orbit contribution to the evolution (or precession) equations for the spins, the conserved 
integrals of the motion associated with the Poincare invariance of the equations (energy, 
angular momentum and center-of-mass integral), and the near- zone metric regularized at 
the location of the particles. We also obtained the bulk metric {i.e. at any field point 
in the near-zone) but with less precision in the 00 and Oi components. The results are 
given for general orbits in the frame of the binary's center-of-mass, and reduced for quasi- 
circular orbits, which is the mostly relevant case for realistic compact binaries detected by 
gravitational waves. We systematically make use of a precise definition of the spin variables, 
which have a conserved Euclidean norm. 

The metric in the near zone is especially of interest in applications such as the numerical 
study of a possible astrophysical environment surrounding black-hole binaries with spins 
[39I , and the comparison with numerical computations of the gravitational self force acting 
on a particle orbiting a Kerr black hole, which are based on first (or second) order black 



hole perturbation theory [40|, |41|. More generally, all the present analytical results (given 
explicitly as they are) should be useful for the accurate comparisons between PN predictions 
and the full numerical computations of the coalescence of two spinning black holes. 

The next step in this current program will be to obtain the next-to-next-to-leading order 
spin-orbit contributions in the gravitational-wave phasing of spinning compact binaries. The 
GW phasing including high-order spin effects (and especially high-order spin-orbit effects) is 
of crucial importance for the parameter estimation of signals received from compact binaries 
by the networks of ground-based as well as, in the future, space-based GW detectors. 
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Appendix A: Comparison with the ADM formalism 

In this Appendix, we give the 3PN correspondence between our spin variables and those of 



Ref. [38[. The 2PN correspondence had been investigated in [29[. In Section VII D of Paper I, 
we showed that our results for the equations of motion were equivalent to the ones obtained 
there in ADM variables, and we constructed the extension of the contact transformation 
required to relate the two methods, as well as the 2PN link between the ADM spin variables 
and the harmonic-coordinates spin tensor. We shall use here the same notations, with an 
overbar for all ADM quantities as well as the convenient short-cut vf = p/m. Let us now 
extend this correspondence in the dynamics to the 3PN spin precession equations. 
The ADM conserved- norm vector S obeys the precession equation 

^ = nxs, (Ai) 

where the precession vector Q can be read directly on the spin-orbit Hamiltonian H^o = 
fli ■ Si + Q2 ■ S2 in Ref. 2^. The two conserved-norm spins start differing at 2PN order, 
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and from the fact that they must have the same magnitude we may write 

S(x,p) = S + 0(x,p) X S + O 



(A2) 



where 6 represents a small rotation vector starting at 0{l/c^). Its expression is found by 
comparing the precession equations for the two spin variables. Indeed, if we plug this relation 
into the precession equation for S, noticing that 0(x, p) = + 0{l/c^), we obtain 



(A3) 



Since we are looking for an expression of that does not feature any spin variable, we can 
forget about the degeneracy remaining along the direction of S and we write 
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Using the method of undetermined coefficients for the dimensionless vector 6, we obtain a 
unique solution (for particle 1), 



with the following expressions for the 2PN and 3PN contributionsH 
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The existence of such a solution for 6 shows also that the 3PN precession equations are 
equivalent in the two formalisms, which gives us further confidence in all results. 



^ Notice that the 2PN terms arc different from equation (6.4) in 
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variable to the different conserved variable Sggp used in Ref. 



29|, since they compare the ADM spin 
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Appendix B: Conversion from the conserved spin vector to the spin tensor 



In this Section we use the same notations for the orbital variables as in Paper I, except for 
the spin tensor, denoted here S^^ . From Eqs. (12 ■2p and f l2.15p . we obtain for our conserved 
norm spin vector, expressed in terms of the spin tensor 5**-', the following expressions: 
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As stated above, these conserved spin variables turn out to be the same as the conserved 
spin vectors of 4^, at 2PN order, when evaluated in the center-of-mass frame. 



Appendix C: Center-of-mass position 

In a general frame, the expressions of the whole set of conserved integrals of motion are 
too lengthy and we therefore chose to present them directly reduced to the center-of-mass 
frame. However, in this Appendix, we provide the complete expression for the center-of-mass 
position G (necessarily in a general frame) since this might be used in the future to reduce 
higher order results to the center of mass. Denoting the PN expansion of G in the form 
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where we have defined for convenience y^ = YaI'^xi- Let us emphasize again that these 
two contributions are the only ones required to derive Eq. (13.31) . i.e. to perform our center 
of mass reduction. In addition to these, the next-to-next-to-leading correction is given by 
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This result translates into a next-to-next-to-leading order correction to the vector z in 
Eqs. fl3.3p . of the form 
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To conclude this Appendix, let us show explicitly that the above (9(l/c'') correction to 
yi written in terms of the CM variables cancels out when reducing our general frame results 
to the center-of-mass frame. Let us denote by z the (9(l/c^) correction to the Newtonian 
expressions of yi and y2, as defined by Eqs. f l3.3p . Since we want to obtain expressions up 
to (9(l/c^) in the CM frame, the terms of order (9(l/c^) in z can only play a role in the 
reduction of the Newtonian parts. These Newtonian contributions are very simple, so we 
can check that the expression of z is needed only up to (9(l/c^). 

The simplest case is that of the relative Newtonian acceleration = — G'mni2/r^2; 
which receives no higher order correction when reexpressing it in the CM frame. The same 
argument applies to the near zone metric. By contrast, the Newtonian energy 



1 2 ( dz\ Gmv 
En = -muv + m — 

2 \dt J 2r 



(C3) 



receives a PN correction when expressed in terms of the CM quantities. However, since the 
terms linear in z cancel out while dz/dt = 0{l/c^), we only need z up to the order 0{l/c^) 
to get the energy up to the 3.5PN level. A similar argument applies to the Newtonian 
angular momentum, which is deprived of linear-in-z contributions: 



rnvx X V + mz X 



dz 

dt 



(C4) 



Appendix D: Orbital angular momentum 

The conserved total angular momentum (such that dJ/dt = 0) is conventionaly split into 
an orbital angular momentum, which will also receive contributions from the spins, and the 
sum of spins of the system: 



L + 



(Dl) 



This actually represents the definition of the orbital angular momentum L, which depends 
on the choice of variables for the individual spins. With this definition, the structure of the 
orbital part is 
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(D2) 



The non-spin terms up to 3PN order can be found in Ref. [5l|- Again neglecting the spin-spin 
terms, we arrive at the following spin-orbit contributions. At leading and next-to-leading 
order, we have 
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As in the case for the energy, the next-to-leading term does not have the same expression as 
the one published for instance in Refs. 23|, |28| , due to our different choice of spin variables; 
but now the difference also shows up in the leading-order spin contributions to L since it is 
actually the next-to- leading order correction in the total angular momentum J. Finally, at 
next-to-next-to-leading order we have 

£3.5PN = (^-^^4 + (^-^^^1^ + ^'-"^^tS^- + ^'-'^^3^ , (D5a) 



35 



= (StO V— (- 
m \ 



19 75 



16 8 



-V 



16 



ij'']v^ + {Sv) V 



+ S 

(3-5)^^ = (^S) n 

+ 



35 65 



z/ + — z/^ h'" + s 



16 16 16 
5m /15 
m 



5m / 5 19 



5 35 65 . 

\ V 

16 16 16 / 



75 



m V 16 8 



V H V V 



16 



75 



z/2 {nvY + 



5m / 15 



m 



5m fi 21 77 2 \ 4 
V 1 ir 

m 



+ (St;) n 



5m ( 51 15 



, — V H V 

m \16 8 



5m / 3 23 



m 



243 



H I — o + TF^ ) (^'^)'^ 



16 



+ (nS*) n 



45 135 



— z/ —V \ (nv) 



33 33 



4 



+ (^w) n 
+ (nS) V 
+ (S?;)v 
+ (nS) V 
+ (5tO V 
+ S 



4 

3 9 



— u H z/ (nv) V + I u u V 



33 83 



11 



33 



4 

5m ( 21 



- + -z/ + — z/^ (nw)'' + — - + —V + — z/ 



17 



109 



8 16 



\nv)v 



m 
5m 



16 



m 



93 



19 











\ [nvf + — 








/ m 




16 y 





5m ( 21 59 



— V (nvy H \ V V 



m. 



35 



15 



-v + 27z/" (nv 



17 295 75 , 

1 z/ 

4 16 4 



43 



{nvY + 



21 33 



+ S 



V 



5m ( 33 



3 31 

V H — V I i;" 

8 8 

9 21 



) {nvf + ( 3 - -z/ - — z/^ ) {nvyv^ + ( "o + — 



9 47 



m 



9 



z/ + -v^ {nvy + 



5m /39 



m 



-z/ - 17z^^ ) {nvfv'^ 



+ ■ 



5m /21 33 



47 

-V + — I V 



(3-5)£2 = (nS) n 
+ (nS*) n 



m V 8 4 

5m / 469 51 



m 



-z/ z/ 

4 



5m / 7 149 



\nv] 



33 



973 39 



m \4 
55 311 



-V H z^^ 

4 4 ' 



V - —u^ invY + —V + Qu^ \v 



4 



5m / 17 103 



m 
5m 



25 



+ (St>) n I — — H — —u — — z/ ) [nv] + (S't;) n I — — H — — z/ — — z/ 



29 247 31 



4 



m 



177 



+ (nS) V -- + —u + 3z/^ (nw) + (Sw) v 



5m f 1 257 



m 



21 2 
-z/ + — z/^ 
4 



+ {nS) V 



117 731 



27 



+ — — z/ H z/^ (nt;) + (S*!;) v — 



25 571 



4 



-z/ H — z/ 
2 



+ S 



17 163 



35 



T^V + —V' 



(nv) 



41 373 



13 



4 



-z/ V 



36 



+ s 



5m 
m 



43 



33 



5m 



m 



-V 



2u 



+ S- 



.5m 



m 



2u 



\nv] 



5m 
m 



153 



31 



-V 



{nS)n 



13 



-u - 2u' 



3 
2 



13 



2u' 



(D5b) 



Appendix E: Gauging away the 3PN spin-orbit terms 

The metric regularized at the location of the two particles, as was presented in Sec- 
tion UllEl displays terms which are formally of 3PN order, that is to say of order 0{l/c^) 
in goQ and 0(l/c^) in g^i. These terms are given by the expressions fl3.11cp and f l3.12cp : 

(for)^=^{(-,S,.) (2-10^) M-24M(-,5,^0}, (El) 
\ =^|2(^,)(nxS)^-^(vxS)4. (E2) 



s"" J ^ c^r^ 



In this Appendix, we show that they are actually a pure coordinate effect, in the sense 
that they can be eliminated by a gauge transformation. We already encountered such a 
gauge transformation in Ref. [49| and in Paper I (see the footnote in Section VI B there), 
where it was used to remove the C(l/c^) terms appearing in the acceleration. Here we 
generalize this transformation to cancel as well the 3PN terms in the components of the 
regularized metric. Notice that, for qqq and goi, we do not control the metric components 
in the bulk at the corresponding order. The transformation we present here is only suitable 
for removing these terms in the regularized metric. 

We consider a coordinate transformation x''^ = + C,^, which has to satisfy = (at 
dominant order) to respect the harmonic gauge condition. We will need to have = (9(l/c^) 
and ^* = 0{l/c^). Since the components gij of the metric must remain unaffected in the 
bulk at order 0{l/c^) (see (13.161) ). we further require that the component C,^ is a mere 
function of time, i.e. does not depend on the field point x. We allow a spatial dependence 
in the component but the harmonicity condition A,^° = 0{l/c^) constrains its structure. 
Hence, we take the following Ansatz for the transformation: 
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A + B^- + B2- + C\d, (-\ + C!,d, (- 
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(E3a) 



e = , (E3b) 



C6 



where A, Bi, B2, CI, and are all mere functions of time. This gauge transformation 
has to be symmetric by exchange of the two bodies, which relates Bi to B2, CI to C^ and 
constrains the expression of A and D^. We could have considered higher order derivatives 
of 1/r in but we found that the above structure was sufficient for our purpose. 

With our simple choice for P, t he gauge transformation changes the metric components 
according to, (see Eqs (6.9) in |46[): 

^€^700 = 2doe + 0(l/c^°) , (E4a) 
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Sigo^ = -doC + ^^e + 0{l/c^) , (E4b) 
S^g.j = -<9.e - d,e + 0(l/c') , (E4c) 

while the acceleration transforms as (see Eqs. (4.8) in [Hij for more general expressions): 

d2 



(E5) 



This shows that ^* is the same as we used in Paper I, and which was called SX^ there. 
Solving for the other unknown functions of time by applying (]E4p to the structure (]E3p . we 
find the following expressions: 
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o 3 
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(E6a) 

(E6b) 
(E6c) 
(E6d) 



This solution allows us, within the class of harmonic coordinates, to cancel the 3PN terms 
(lEip in the regularized metric, which are therefore pure gauge. 
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